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Abstract 

We investigate theoretically the energy levels of two identical Fermions in a harmonic potential 
well under p-wave Feshbach resonance. The magnetic fields needed to affect the energy levels differ 
strongly among the levels, and they can be at values far from the free space resonant field. 
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There are substantial recent activities in studying cold atomic gases in the presence of 
an optical lattice (e.g. [l]). Part of this interest stems from the possible connection with 
solid-state physics (e.g. |2]) and the vast theoretical literature on exotic states in lattice 
hamiltonians. Generally, the interaction between atoms are weak. To increase the interaction 
between particles, one possible way is to make use of Feshbach resonance. Indeed, a few 
such experiments have already been carried out 3|, Ij]. An important question arises already 
even when the tunneling of atoms between different potential wells can be neglected. For 
more than one particle in a well, the energy levels, wavefunctions etc are strongly modified 
by the Feshbach resonance. Theoretical studies of this "dressed state" have been carried 
out in some recent papers j^, ^. Experimentally, Ref ^ has been able to obtain some 
important information about this dressed state, namely, the probability of finding a particle 
in a particular harmonic oscillator level. They did this by starting with a deep lattice 
potential and, by sweeping the magnetic field through resonance and then turning off the 
optical lattice, occupation of the different Bloch bands was experimentally observed. The 
energy of the lowest bound state has also been measured by radio frequency spectroscopy 
4]. Understanding these dressed state is also crucial for obtaining the lattice Hamiltonian 
3], the first essential step to eventually able to understand any many-body systems to be 
studied. 

The above references are on s-wave Feshbach resonance between two Fermions 

of different internal hyperfine spins. Here we would like to extend the study to two iden- 
tical fermions (same hyperfine spin) under a p-wave Feshbach resonance. P-wave Feshbach 



resonances in Fermions (e.g. [8|, l9|, llO| ) itself have generated much activities. One reason is 
possible p-wave superfiuidity (e.g. [ll, 12,[l^, Q). There has also been a suggestion 15] to 



obtain molecular superfiuid by first starting with a lattice with two atoms per site and then 
"melt" the Mott state by turning off the optical lattice. This may perhaps be a route to 
obtain p-wave superfluidity and gives additional impetus to study Fermionic atoms under 
p-wave Feshbach resonance within an optical lattice. 

In this paper, we shall then consider two identical Fermions in a harmonic potential 
well under a p-wave Feshbach resonance. We shall flnd how the energy levels depend on 
the magnetic fleld detuning. We shall show that, in contrast to s-wave resonances, due to 
the strong momentum dependence of the effective interaction, there is a large shift in the 
magnetic fleld from resonance which is necessary to affect the energy levels with the shift 
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itself depending strongly on the energy level under study. The validity of this prediction 
can be tested by performing an experiment similar to that of 3j. We employed a method 
different from the existing literature, an approach some readers may find more transparent. 
We shall also report on the generalization to resonance with higher angular momenta. 

In a harmonic trap, we can first separate the center of mass versus the relative motion, 
and we need only to treat the latter. The wavefunction of the "dressed" molecule in the 
relative coordinate r consists of two parts, ilj{r) in the open channel and '-)ipc{r) in the closed 
channel. Here, for later convenience, we have introduced a normalized wavefunction v'c(^ 
in the closed channel and denote the amplitude of our dressed molecule in this channel by 
7. The normalization of the full wavefunction is thus I7P + = 1, where is the 

short-hand for j d^r. For p-wave resonance, due to the dipole interaction, the resonant 
field depends on the angular momentum mi along the applied magnetic field direction. For 
example, for ^'^K in the ||,— | > hyperfine state, the = resonance occurs at a field 
around 0.5G above the mi = ±1 ones |^. For definiteness, in below we shall assume that 
we have the mi = resonance only. Modifications need to other m^'s are evident. 

The Schrodinger equations are 

+ U{f) j ij{f) + W{r)-fyDc{r) = Ei){r) (1) 
W{r)^{r) + v'^ipc{r) = E^^,{f) (2) 

Here = m/2 is the reduced mass, U{r) is the trap potential, and ^pd'r) is the closed- 
channel molecule of / = 1, m; = symmetry with (unrenormalized) energy v\ and W{r) 
is the coupling between the open and closed channels, and E is the energy of the dressed 
state. Here, we have already truncated the Hilbert space for the closed-channel to one 
single state responsible for the resonance of interest. For simplicity, we have also suppressed 
the part of the wavefunction representing the internal (hyperfine) degrees of freedom which 
distinguishes between the open and closed channels. Eq ([I]l2]) are the same as those written 
down for the s-wave (e.g. p!]) except the symmetry of V'c(^- 

Equation ([1]) can be solved in terms of the Green's function defined by 

f — + f/(f) - e\ GE{r, ^) = —5%r- f") . (3) 
Y 2mr J mr 

Substituting the result into eq ([2]), we obtain the implicit equation for the energy E as 

E-u' + — [ ^lif^GEif, n^f") = (4) 



where we have introduced the short-hand <^c(^ = W{r)(pc{''^- The problem thus reduces to 
finding GE{'r,r'). We here adopt a different method from ref 16|. We would also like to 
express our final answer in terms of experimentally available quantities. For the latter, we 
first observe that, in the absence of the trapping potential, the expression in eq (jl]) is related 
to a quantity that appears in the scattering amplitude. Explicitly, the scattering amplitude 
for incoming (outgoing) wavevector k {k') can be found to be 



(5) 



ik-r 



where G^^^'^{f,f') obeys eq ([3]) with U{r) set to zero. Here E = ^pc{k) = ^^(pc{f)e' 
For scattering in the Z = 1, = channel, f^{k') = 4:TcfiQ(k)Y^*{k)Y^{k') where /lo has 
the following form at small energies: 



/lo(fc) 



(6) 



- - + ck"^ — ik^ 

V 

Here v has the dimension of a volume and c inverse length. These parameters, in principle 
available from experiment (e.g. [s!]) or other theoretical calculations, will be used as an 
input below. Near the resonant field Sq, — ^ is roughly linear in {B — Bq), whereas c can 
be regarded as a constant. 

At small energies, <^c(^) ~ —ik^J.^^Pci^^z = —ik^y/Sa, where we have also defined the 
coupling constant a (This definition is identical with that in jisl . We then have the 



correspondence 



rrir 



free. 
E > 



f,r')ipc{r') 



- + ck^ 

V 



ik^ 



(7) 



This equation, when expanded and comparing powers of fc, yields the same "renormalization" 
relations as those written before in momentum space in, e.g. 
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The Green's function GE{r,r') and G^^^'^{f^r') are both singular as r ^ r*, diverging as 
pr^, but otherwise finite. Defining G^^^{f,r') = GE{r,r') — G^^^'^ {f, r') , G^^^{r,r') is thus 



finite for all r and r*. The equation ([ 

-- + ck^- ik^ 

V 

It can be verified that obeys 



for the energy can then be written as 



\a\ 



2m, 



V + f/(f) -E]G 



■^reg 



r, r 



-U{r)Gi 



r, r 



(8) 



(9) 



Below, we shall find G'^^^{r,r') for the case of a spherical trap of frequency cu, thus 

U{r) = ^nirUj'^r'^- It is sufficient to do this for r > r' since the r < r' part 

can be obtained by the symmetry of G^^^{r,r'). We have G^j^^^{f,r') = ' = 

'i'lrik^i ji{kr')h\^\kr)Yirni{r)Yi^^{f') for r > r'. By the symmetry of (fc, we see that 
only the terms with Z = 1, m; = would contribute. Writing this part of G^^^{f,r^) as 
gi{r,r')YiQ{f)ji{kr')Y*Q{f'), we are thus left to solve 

^ ^ r^4- + + f^(^) - ^1 = -f/(r)47riA;/i;')(A;r) (10) 



(ir dr 2mrr'^ 

Fortunately, a particular solution for eq ( fTOl) can easily seen to be —A'n-ikh^^\kr). The 
origin of this is simply that — G';^'^^^(r, r') is a particular solution to eq (l9l). The homogeneous 
solution to eq ffTOl) can be expressed in terms of the confluent hypergeometric functions F, 
so we obtain thus 



gi{r, r' 



Ur\^ . (1) 



e -Anikh\ \kr) 

(11) 

where /ii = ^ — |, /i2 = ^ + | = /Ui + |, and Ir = ^Jhjrr^w is the oscillator length for 
relative motion, and &i(r') and 62('"') are coefficients to be determined. The condition that 
Gsir, r*) — > as r — ^ oo gives bi ^l-^^) + ^27^z^f^ = 0, that is bi = |f(E]^^2- For small r', 
ji{kr') — > kr' /?). The condition that GE{r,r') must be regular in the limit limr'^o hnir-^rY 
therefore requires that there cannot be any term in gi{r, r'). Using h^l\kr) —i/{kr)'^ 
as r 0, we therefore have limr'^o b2{r') = On the other hand, for small r, we obtain 
from eq ffTTl) that gi{r, r') = bi{r-')j- — ^ik^r. (It can be verified that, with the above value 
of 62, there are no terms proportional to or in gi{r,r')). Substituting this expression 
into the right-hand-side (RHS) of eq ([H]), we see that, in the limit of short-ranged (fc, we 
need only the value of bi{r' 0), and the RHS can be evaluated to be - ik^. The 

—ik^ terms on the two sides of eq ([8]) cancel. Using the above obtained value for 61 (0) and 
the definition for /xi and fi2, we obtain finally the equation 

^ I rP - ^ I'd ~ (-,0) 

for the energies of the dressed states, expressed entirely in terms of quantities entering the 
scattering amplitude ([6]) and the harmonic trap frequency u {k^ = 2mrE). 
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It is convenient to rewrite this equation in dimensionless form: 

- ^ - (-2c/,)- = 8 ;\ (13) 

For "^^K ||, — I >, — c ^ 0.02ao \ whereas ~ 1600ao for a trap of frequency u = 70kHz, 
a typical experimental number. Thus the product (— 2c/r) is large and ~ 64. The resulting 
energy levels are plotted in Fig [H For "high" magnetic field, (— ^ —>■ oo), the energy 
levels are given by the non-interacting values (2n + ^)u! {n is a non-negative integer). For 
= 0, the state has energy decreasing with field, becomes linear in —1/v for — < 100, 



and thus becomes indistinguishable with the the free space bound state energy H 



th = — (-cv) ^"^^ B < Bq. For n > 1 and decreasing field, the energy of the level shift to 



(2n + i)u; = (2(n — 1) + 1)^7. However, the value of the field where the energy crosses from 
its high field value to the low field one depends on the level n, increasing as n increases. The 
width of field for this transition also increases with n. This is due to the large (— c/,) value 
discussed above. 

The energy levels can be measured by radio-frequency spectroscopy ^]. Below, we analyze 
instead the implication on an experiment similar to that in [si] (theoretically analyzed in 
0, Q), summarized already briefly in the introduction. We thus evaluate, for a given 
energy level of Fig [T|, the probability of finding a particle in the single particle harmonic 
oscillator states. For this, we first have to express the twcnparticle wavefunction ilj{f) in 
harmonic oscillator basis. A procedure similar to that in joj gives ip{r) = ^nVn'^nioi^)-, 
where 

{lS# + E„[A./(£-(2n + i))]'^} 

'^niQ^r) = 3^"5/2 ^cos6'F(— n, |, 1^) )e~2(i;:) are the normalized wavefunctions with I = 1, 



mi = 0, An 



3/4«^/2 ^^^^-^ V '"^ 2' \l 

1 (2n+3)!! 



3 2"-in! 



ll/2 

are normalization coefficients. The amplitude in the closed 



channel is 7 = [1 — l^nP]^^^- Below, we consider the special case where the term oc m 



the denominator of eq ffT4l) is negligible, correspondingly 7—^0. (This assumption affects 
only quantitatively the probabilities at intermediate fields given below but not the qualitative 
conclusions.) Using this iplr) and assuming that the center of mass motion is in its harmonic 
oscillator ground state, it is straight-forward to express the resulting wavefunction in the 
coordinates fi, r2 of the two particles. For the latter wavefunction (p, we shall employ 
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instead the quantum numbers Ux, Uy and since these are more directly related to the 
occupation of the different Bloch bands and thus the momentum states after turning off of 
the optical lattice. For example, the state with quantum numbers n^,, n^, contributes to 
momentum states with nx[y^z)2 < \kx{y,z)\ < {i^x{y,z) + l)f ; where d is the distance between 
neighboring sites in the lattice [6.]. The resulting expressions however are lengthy and we 
shall not write them down here, but only give the probabilities below. For definiteness, we 
consider an initial field B with initial two-particle state corresponding to relative 

motion with quantum numbers n = 0,/ = l,m = (Energy of relative motion = ^u.) In 
the two-particle basis, the state is given by ^ [0000(^1)0001(^^2) — 0ooo(^2)0ooi(n)] where the 
subscripts stands for {nx,ny,nz) = n. The initial probabilities of finding a particle in the 
single-particle state 0^ are given by Pqoo = -Pooi = 1- ( -Pn = for all other n). The final 
probabilities of finding a particle in state n when the magnetic field is swept to the final value 
B where the corresponding scattering parameter is 1/v is plotted for a few n's in Fig[2l The 
transition from the "low" field values to the "high" field values occur around —l^/v ~ 200 
(c./. Figd]) for this particular state (This field increases with n). For magnetic fields with 
—1/v above this value, the probabilities correspond to the two-particle relative motion in 
the = 1, / = 1, m = state. We find Pqoo = i, Aoi = ^002 = ^, ^100 = ^010 = 2^, 
-Pioi = -foil = -P003 = J), with all other = 0. The field at which the probabilities 
switch from one set of values to the other can be used to indicate when the energy levels 
change from the low to high field values. For in a trap with 70kHz considered above, 
the shift ll/v ~ —200 corresponds to — - ~ 0.5 x lO^^ao'^, or a field ^ 1.15G above the free 

nn 

space resonance [8|, 118|, a large shift that should be readily discernible by experiments. This 
shift in field is roughly proportional to 1//^ hence uj. 

Finally, we discuss Feshbach resonances in higher angular momentum channels, say I, 
mi (Bosons or non-identical Fermions if / even). Eq ([5]) is still valid, and the scattering 
amplitude has the form (generalizing eq ([6])) A;^'/ ( — ar+r + Qm^^ + ... — i/c^'^^ j , where aim 
has dimension of a length. One have relations similar to eq (JTj) and (IHl) with coupling 
constants aim = [■^^yA Jf'fcir)r^Yi^{f). We can obtain the corresponding part of G^^^ with 
relevant /, mi symmetry with gi{r,r') given by an expression similar to eq ffTTj) except that 
the bi{r') (&2(^')) term has the factors ( ( f- ) ) with appropriate changes in the 

arguments of the confiuent hypergeometric functions F, and h^i^ . If we proceed as 
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in text and include only the regular term oc r' in gi{r,r'), we would get 

-, / r, N 21+1 p/2^h3 _ E_\ 

-4n + --^' + - = (-i)'"Mf w^l-i (15) 

«/m VrJ T{ — 

This formula reduces to the known / = results in the literature 



161] and the / = 1 

result above (eq (|T2l) ). However, this is incomplete if / > 2. For / = 2 for example, the 
term oc 62? together with the in2{kr) term from /;,2^''(A;r), produce a term oc ^ x = r 
as the lowest order contribution to g2{f,r')- (It can be verified that there are no terms of 
lower order in r). As a result, the RHS of eq (ITSl) acquires an extra term i?2m given by 



2r(§)/^^,(r>y,*^(f) 



«2n 



The integral involves short distance behavior of ipc, and in contrast to 



the coupling constant aim, cannot be eliminated from the final answer using the knowledge of 
the low energy scattering amplitude. Similar "residual" terms Rim arise for all / > 2, though 
the expressions become more involved and will not be given here. However, fortunately, for 
a given Im, Rim is independent of the energy, magnetic field (if near resonance) and the 
energy level under consideration. Since near the resonance, ^j+r is expected to be linear in 
the deviation of B from resonance, the extra term Rim simply gives a shift of the resonant 
field common to all the energy levels for given Im and trap parameters Ir- Apart from 
this, the behavior of the energy levels as a function of field is thus rather similar to that 
depicted in Fig [T] (though with energy levels given by {2n + ^^)c<j far from resonance). The 
quantitative details will of course depend on the values of the parameters cim etc. 

In conclusion, we have calculated energy levels of two Fermions in a harmonic potential 
well under a finite angular momentum Feshbach resonance. In contrast to wide s-wave 
Feshbach resonances, there are large shifts in the resonant fields dependent on the energy 
evels. This prediction can be tested by experiments similar to those already performed in 
3,3. 
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FIG. 2: The probabilities P of finding a particle in the state n for n = (000) (solid), (001) (dashed), 
(100) (long dashed), and (002) (dot-dashed), under the field sweep procedure described in text, v 
is unit of /j?. 
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